The Effects of Next-Nearest-Neighbor Interactions on the Orientation Dependence of 
Step Stiffness: Reconcihng Theory with Experiment for Cu(OOl) 
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Within the solid-on-solid (SOS) approximation, we carry out a calculation of the orientational 
dependence of the step stiffness on a square lattice with nearest and next-nearest neighbor inter- 
actions. At low temperature our result reduces to a simple, transparent expression. The effect of 
the strongest trio (three-site, non pairwise) interaction can easily be incorporated by modifying the 
interpretation of the two pairwise energies. The work is motivated by a calculation based on nearest 
neighbors that underestimates the stiffness by a factor of 4 in directions away from close-packed 
directions, and a subsequent estimate of the stiffness in the two high-symmetry directions alone that 
suggested that inclusion of next-nearest-neighbor attractions could fully explain the discrepancy. As 
in these earlier papers, the discussion focuses on Cu(OOl). 



PACS numbers: 68.35.Md 81.10.Aj 05.70.Np 65.80.-l-n 

I. INTRODUCTION 

At the nano-scale, steps play a crucial role in the 
dynamics of surfaces. Understanding step behavior is 
therefore essential before nano-structures can be self- 
assembled and controlled. In turn, step stiffness plays a 
central role in our understanding of how steps respond to 
fluctuations and driving forces. It is one of the three pa- 
rameters of the step-continuum model^ which has proved 
a powerful way to describe step behavior on a coarse- 
grained level, without recourse to a myriad of micro- 
scopic energies and rates. As the inertial term, stiff- 
ness determines how a step responds to interactions with 
other steps, to atomistic mass-transport processes, and 
to external driving forces. Accordingly, a thorough un- 
derstanding of stiffness and its consequences is crucial. 

The step stiffness (3 weights deviations from straight- 
ness in the step Hamiltonian. Thus, it varies inversely 
with the step diffusivity, which measures the degree of 
wandering of a step perpendicular to its mean direction. 
This diffusivity can be readily written down in terms of 
the energies Sk of kinks along steps with a mean ori- 
entation along close-packed directions ((110) for an fee 
(001) surface): in this case, all kinks are thermally ex- 
cited. Conversely, experimental measurements of the 
low-temperature diffusivity (via the scale factor of the 
spatial correlation function) can be used to deduce the 
kink energy. A more subtle question is how this stiffness 
depends on the azimuthal misorientation angle, conven- 
tionally called 9 and measured from the close-packed di- 
rection. In contrast to 9 = steps, even for temperatures 



much below Sk, there are always a non- vanishing number 
of kinks, the density of which are fixed by geometry (and 
so are proportional to tan 9). In a bond-counting model, 
the energetic portion of the step free energy per length 
(or, equivalently, the line tension, since the surface is 
maintained at constant [zero] charge^) f3{9) is cancelled 
by its second derivative with respect to 9, so that the 
stiffness is due to the entropy contribution alone. Away 
from close-packed directions, this entropy can be deter- 
mined by simple combinatoric factors at low temperature 

3.4.5 

Interest in this whole issue has been piqued by the re- 
cent finding by Dieluweit et al& that the stiffness as pre- 
dicted in the above fashion, assuming that only nearest- 
neighbor (NN) interactions ei are important, underesti- 
mates the values for Cu(OOl) derived from two indepen- 
dent types of experiments: direct measurement of the 
diffusivity on vicinal Cu surfaces with various tilts and 
examination of the shape of (single-layer) islands. The 
agreement of the two types of measurements assures that 
the underestimate is not an anomaly due to step-step 
interactions. In that work, the effect of next-nearest- 
neighbor (NNN) interactions €2 was crudely estimated 
by examining a general formula obtained by Akutsu and 
Akutsuji showing a correction of order exp{—e2/kBT), 
which was glibly deemed to be insignificant. In subse- 
quent work the Twente groupS considered steps in just 
the two principal directions and showed that if one in- 
cluded an attractive NNN interaction, one could evalu- 
ate the step free energies and obtain a ratio consistent 
with the experimental results in Ref. 1^. This group later 
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extended their calculations^^ to examine the stiffness. 

To make contact with experiment, one typically first 
gauges the diffusivity along a close-packed direction and 
from it extracts the ratio of the elementary kink energy 
Ek to T . Arguably the least ambiguous way to relate Sk 
to bonds in a lattice gas model is to extract an atom from 
the edge and place it alongside the step well away from 
the new unit indentation, thereby creating four kinks. 
The removal of the step atom costs energy 3ei-|-2e2 while 
its replacement next to the step recoups ei + 2e2. Thus, 
whether or not there are NNN interactions, we identify 
£fc = — ^ei =5|ei| (since the formation of Cu islands im- 
plies ei < 0); thus, as necessary, Sk > 0. Note that for clar- 
ity we reserve the character e for lattice-gas energiesji^ 
which are deduced by fitting this model to energies which 
can be measured, such as Sk- 

The goal of this paper is to compute the step line ten- 
sion (3 and the stiffness /? as functions of azimuthal mis- 
orientation 9, when NNN (in addition to NN) interac- 
tions contribute. Since it is difficult to generalize the low- 
temperature expansion of the Ising model, we instead 
study the SOS (solid-on-solid) model, which behaves very 
similarly at low temperatures and at azimuthal misori- 
entations that are not too large, but can be analyzed 
exactly even with NNN interactions. This derivation is 
described in Section II, with most of the calculational de- 
tails placed in the Appendix. In Section III we derive a 
simple expression for the stiffness in the low-temperature 
limit, presented in Eq. (|14|l . We also make contact with 
parameters relevant to Cu(OOl), for which this limit is 
appropriate. In Section IV we extend the formalism to 
encompass the presumably-strongest trio (3-atom, non- 
pairwise) interaction, showing that its effect can be taken 
into account by shifting the pair energies in the preceding 
work. The final section offers discussion and conclusions. 



II. NNN SOS MODEL ON A SQUARE LATTICE 

Including NNN interactions in the low-temperature ex- 
pansion of the square-lattice Ising model lifts the remark- 
able degeneracy of the model with just NN bonds. In 
that simple case, the energy of a path depends solely on 
the number of NN links, independent of the arrangement 
of kinks along it; thus, the energy of the ground state 
is proportional to the number of NN links of the short- 
est path between two points, and the entropy is related 
to the number of combinations of horizontal and vertical 
links that can connect the points. '^•^ Including NNN in- 
teractions causes the step energy to become a function of 
both the length of the step and the number of its kinks, 
eliminating the simple path-counting result,^ It can then 
become energetically favorable for the step to lengthen 
rather than add another kink. This causes the NN en- 
ergy levels to split in a non-trivial way, making it possible 
for a longer step to have a lower energy than a shorter 
step. A related complication is that the expansion itself 
depends on the relative strength of the NNN-interaction: 



Instead of an expansion just in terms of exp(— jeil/fesT), 
the expansion also is in terms of exp(e2/2/cB2^)- Hence, to 
take the NNN-expansion to the same order of magnitude 
as the NN-expansion, an unspecified number of terms is 
required, depending on the size of the ratio 62/61. 

Since the NNN Ising model cannot be solved exactly 
and we cannot generalize the low-T expansion, we turn to 
an SOS model, which was used in earlier examinations 
of step problems, most notably in the seminal work of 
Burton, Cabrera, and Frank,-!^ and later used for steps of 
arbitrary orientation by Lcamy, Gilmer, and Jackson. ^4 
It was also applied to an interface of arbitrary orientation 
in a square-lattice Ising model. 

Although the SOS model can be treated exactly, 
the result is somewhat unwieldy. Fortunately, at low 
temperature — the appropriate regime for the experi- 
ments under consideration — the solution reduces to a 
simple expression. 



A. Description of Model 

Consider a step edge of projected length L separating 
an upper adatom-free region from a lower adatom-fiUed 
region (see Fig.^. The step edge is completely described 
by specifying its height yi at position i {0 < i < L). The 
energy of the step edge depends on the number of broken 
bonds required to form it. Let V and H represent the 
vertical and horizontal NN bond strengths divided by 
kgT, and let U and D represent up-diagonal and down- 
diagonal NNN bond strengths over fc^T. Then the step- 
edge energy E = E{{Ai}) depends only on = yi — 
Vi-i- 

For clarity, we consider two examples. First, if A; = 3 
(as is the case between columns a and b in Fig. then 
between positions i and i + 1 there are 3 broken i7-links, 
2 broken [/-links, and 4 broken D- links. There are also 
2 broken T^-links, but this number is independent of A^, 
since every step-edge configuration of projected length L 




FIG. 1: A finite-sized step edge whose projected lengtli is L. 
The step has height yi at position i {0 < i < L). The height 
difference ul — yo is fixed; thus, the step edge makes an angle 
9 with the horizontal axis, and has an overall slope m (shown 
as the top of the gray region). The energy of the step edge 
is found by counting the number of broken links required to 
form it. Here all NN and NNN broken links are shown. 
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requires exactly L broken V^-links. Similarly, if = — 3 
(as is the case between columns c and d in Fig. then 
there would be the same number of broken if-links, but 
there would now be 4 broken [/-links and 2 broken D- 
links (that is, the number of broken U and D links switch 
from the previous case). From these examples we see 
that, in general, there are |Ai| broken iJ-links, |Ai — 1| 
broken [/-links, and |Ai -I- 1| broken D-links. It therefore 
follows that the step-edge energy is 



Y,{V + H\^,\ + U\^,-1\+D\X + 1 

^ 2—1 

L 



Because we seek the orientation dependence of /3 and 
/?, we constrain the step to have an overall offset Y = 
Vl — yo = Lta.n9 — X^i^i (This constraint is repre- 
sented in Fig. n by the shaded gray area. Equivalently, 
we specify that the overall slope of the step is to = tan 9.) 
The constrained partition hmction is therefore 



{A} 



,-£;({A.})/fcsT 



(2) 



where {A} is the set of all A^ each of which ranges over 
all integers. From Z{Y) we can find the orientation de- 
pendence of the free energy F{Y) = —kBTln Z{Y), the 
projectedivee energy per length f{m) — F(Y)/L, and the 
line tension (or free energy per length) f3{9) = f{m) cos9 
(since the step length is L/ cos9); thence, we can find the 
stiffness (3(9) = (3{9) + d"^ (3{9) / 89^ . 

For future reference, note that the process of extract- 
ing an atom from the step-edge and replacing it alongside 
the edge, discussed in the penultimate paragraph of the 
introduction, creates two pairs of A = -1-1 and A = — 1, 
costing AH according to Eq. and removing a net of 
2 NN bonds, so that H = -ei/2fcBT = ek/ksT. Simi- 
larly, we compare the energies of two NN atoms, abut- 
ting [the lower side of] a step edge ({A^} = 0) at Iq and 
either parallel or perpendicular to the edge. In the first 
case, Ai(, =-|-l and AiQ+2 = — 1, with an added energy of 
2H + 2{U -\-D) according to Eq. In the perpendic- 
ular case Ai(, = +2 and ^io+i — —2, implying an added 
energy of 47? -|- 4(11 + D). Counting bonds we see that 
the parallel configuration has one more ei bond and two 
more £2 bonds than the perpendicular configuration. In- 
voking H = -ei/2kBT, we see that U + D = -e2/fcBT; 
if [/ = Z?, then D = -e2/2kBT. The factor-of-2 differ- 
ence between broken links in Eq. ^ and broken bonds 
was noted (for H links) already in the classic exposition 
by Leamy et alii^ An alternate argument, presented over 
a decade agOfi^ for this factor of 2 is that the ragged 
edge is created by severing bonds along the selected path 
through an infinite square. This leads to the formation 
of two complementary irregular boundary layers (with 



opposite values of {A^}, so that the associated energy 
of each is half that of the broken bonds (at least when 
U = D). 



B. Evaluation of the Free Energy 

As detailed in the first part of the Appendix, the sum 
in the Fourier transform of Z{Y), which we denote by 
W{ii), factorizes. Thus, it can be written as 

W{^l) = exp [~Lg{ifi)/kBT] , 

where g{i^J.) is the reduced Gibbs free energy per column. 
To evaluate the inverse transform, we exploit the saddle 
point method and obtain (see Appendix for details) 



Z{Y) « exp 



—L pq tan( 



.9(^0) 
knT 



(3) 



where the saddle point (/io = 
by the stationarity condition 



is defined implicitly 



g'jpo) 

knT 



= 171 = tan0. 



(4) 



Here, prime (as in g') denotes a derivative with respect 
to p. This result can be regarded as applying a "torque" 
to the step to produce a rotation 9 = tan~^ m from the 
minimum-energy, close-packed orientation^ 

Taking the logarithm of Eq. Q, we find the projected 
free energy per column /(to) as a Legendre transform of 
the reduced Gibbs free energy per column g{po)- 



knT 



Pom 



gjpo) 

knT 



(5) 



Note that this expression is valid only for L ^ 1; for 
finite-sized systems, corrections are required. As stan- 
dard for Legendre transforms ji^ we have 



knT 



/M = 

ksT g"{po)' 



(6) 



where / = d^f/dm^. Using j3{9)a = f{m)cos9 and 
TO = tan0, with a the lattice constant of the square (i.e., 
the column spacing, which is (l/-\/2) the conventional fee 
lattice constant), we can rewrite the stiffness as 



(3{9)a = /(to)/cos3 6», 
or, similar to results by Bartelt et al.^ 
ksT g"{po 



I3{9)a 



kBT 



■ cos^ 9 



(7) 



(8) 



Thus, we only need g"{p) to find the stiffness as a func- 
tion of to or 9. 

Of course, po in 9" must be eliminated in favor of to via 
Eq. The details for the general case are somewhat 
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involved. Here, we simplify to the physically relevant 
case oiU ^ D and, defining S = H + U + D ^ H + 2D, 
just quote the results: 



iPo) 



2 sinh po 
C(^,Po) 



coth po 



(9) 



where C{S, po) = cosh S — coshpo and po(m) is found by 
inverting 



sinh po sinh S' 



C{S, Po) [sinh 5 - C(5, po) (1 - e-^^)] 



(10) 



Some details can be found in the Appendix. Since 
Eq. Hll)|l is a quartic equation for coshpo or e''°, the ex- 
plicit expression for pa{ra) is rather opaque. However, at 
low-temperatures, a simpler formula emerges, as shown 
in the next section. 



III. LOW-T SOLUTION: SIMPLE EXPRESSION 

At low temperatures, we find that the appropriate root 
for pq diverges. Then we can write coshpo ~ sinhpo ~ 
ePo/2. Of course, H cxl/T so that coshS* « /2. With 
these approximations, Eq. (|1U|I becomes quadratic in e''": 



oPO+S 



m 



(eS _ ePo)[eS - (e^^ - eP")(l - e-^^)] 
Likewise, the expression for g"{po), Eq. becomes 



(11) 



9" (po) 



2ePo 



(12) 



Solving for e^" in Eq. Hll|l and inserting the solution into 
Eq. gives 



9" jpo) 
kuT 



= —fnd (1 — m 



4me 



-2D 



(13) 



so that, from Eq. ((HJ, and recalling D= —e2/2kBT, we 
arrive at our main result, a simple, algebraic expression 
for /3 as a function of m: 



keT 
Pa 



)2 -I- 4me<^2/fcBT 



(1 +m2)3/2 



(14) 



We examine Eq. H14|l in several different limiting cases. 
When £2 = 0, this reduces to 



knT 



m + m 



/3a (l + m2) 



(15) 



as found in a previous study involving only NN 
interactionsi^ Interestingly, at 6* = 45°, Eq. H14|) shows 
a simple dependence on €2, namely. 



keT 
/3a 



(16) 



Of course, this reduces to the venerable Ising result of 
I/V2 in the absence of NNN interactions {e2 = 0)M^ 
By considering just the lowest and second lowest en- 
ergy configurations,^'^'' Zandvliet et al. obtained the 
resul^lfi (expressed with our sign convention for 62) for 
the maximally misoriented case m = 1 



keT 

(3a 



V2 



(17) 



which has, for the attractive £2 of primary concern here, 
some qualitative similarities to Eq. ll()l) (including the 
value 1 / V2 for £2 = 0) but is too small by a factor of 2 for 
t2/2kBT ^ 0; even the coefficient of the first-order term 
in an expansion in £2/2fcsr is half the correct value. For 
the opposite limit of repulsive £2 , Eq. H17|) levels off (at 
\/2), in qualitative disagreement with the actual expo- 
nential increase seen in Eq. H16() . 

Fig. [3 compares Eq. (|14|l to corresponding exact so- 
lutions [found by numerically solving Eqs. lO, ®: and 
(|10|l ] at several temperatures when £2 = £i/10. We see 
that Eq. (|14|l overlaps the exact solution at temperatures 
as high as Tc/6. As the temperature increases, the stiff- 
ness becomes more isotropic, and Eq. (I14fl begins to over- 
estimate the stiffness near 9 = 0°. 

Finally, in Fig. O (using the experimental value^^ Sk = 
128 meV ^ £1 = -256 meV), we compare Eq. (O to 
the NN-Ising model at rj= 320K, as well as to the ex- 
perimental results of Ref. 6. For strongly attractive (neg- 
ative) £2, ksT/Pa decreases significantly. In fact, when 
£2/£i is 1/6, so that -e2/2kBT = {€2/ ei){ekJkBT) w 
(1/6)4.64, the model-predicted value of kgT/Pa has de- 
creased to less than half its £2 = value (viz. by a factor 
of 0.46, vs. 0.63 if Eq. (UTj) is used), so about 3/2 the ex- 
perimental ratio. If £2/£i increases even further, ksT / [3a 
further decreases and develops positive curvature, caus- 
ing an cndpoint local minimum to appear at = 45°. We 
can determine when this occurs by expanding Eq. (|14|) 



keT 



£2 =ei/10 






T = Tc/3 




T = Tc/6 


■ 


T = Tc/10 




low-T = Tc/3 




low-T = Tc/6 




low-T = Tc/10 



m = tan 



FIG. 2: The range of validity of Eq. I|14|l is examined by 
comparing it to exact numerical solutions of the SOS model 
at several temperatures. In the legend Tc refers to the NN 
lattice-gas (Ising) model; for [eij = 256meV, Tc = 1685K. 
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about m = 1: 

Setting the coefficient of (m — 1)^ to zero gives —2D = 
— — ln(6) —1.8, wfiich corresponds to a value 
of ksT/Pa = \/3/6 w 0.29, about 2/5 tfie value at 
£2 = 0. For r = 320K and Sk = 128 meV, this corre- 
sponds to £2/61 « 0.2. However, for the NNN interaction 
alone to account for the factor-of-4 discrepancy between 
model/theory and experiment reported by Dieluweit et 
aL^, Fig. 131 shows that £2/61 ~ 0.3 would be required. 

IV. EFFECT OF TRIO INTERACTIONS 

In addition to the NNN interaction, trio (3-atom, non- 
pairwise) interactions may well influence the stiffness. 
The strongest such interaction is most likely that asso- 
ciated with 3 atoms forming a right isosceles triangle, 
whose sides are at NN distance and hypotenuse at NNN 
separation. In a lattice gas model, there is a new term 
with CRT times the occupation numbers of the 3 sites^ 
Note that this trio interaction energy e^^T is in addition 
to the contribution 2ei-f £2 of the constituent pair inter- 
actions. If wc count broken trios and weight each by i?, 
we find an additional contribution to Eq. 1^ of R times 

4|A,| -f 2(5a..o + 2 = 2|A,| + |A, + 1| + |A,-1| + 2, (19) 

where we have converted the Kronecker delta at i = 
to make better contact with Eq. Thus, without fur- 
ther calculation we can include the effect of this trio by 

keT 



^<»)a T = 320K 




0.2 0.4 0.6 0.8 1 



FIG. 3: Eq. H14^ is plotted for a variety of different val- 
ues of D = -e2/2fcsr, where ei and 62 are NN- and NNN- 
interaction energies, respectively, in a lattice-gas picture. The 
solid curve denoted "Ising NN" corresponds to £2 = 0. The 
dots labeled "Exp't" are taken from Fig. 2 of Ref. (3 and were 
derived from the equilibrium shape of islands on Cu(OOl) 
at 302K, with the line segments serving as guides for the 
eye. To minimize clutter, we omit similar data derived from 
correlation functions of vicinal surfaces at various tempera- 
tures. Note that for £2 =ei/4 a maximum has developed near 
tan 6 = 1/2 that is not evident in the experimental data. 



replacing H hy H + 2R, U hy U + R, D hy D + R, and 
(trivially) VhyV+2. 

By arguments used at the end of Section IIA, we recog- 
nize R — —^eRT- Consequently, the effective NN lattice- 
gas energy is £1 -I- 2ei^T and, more significantly the ef- 
fective NNN interaction energy is £2 + £flT- Thus, cut 
must be attractive (negative) if it is to help account for 
the discrepancy in Fig. 2 of Ref. between the model 
and experiment. Furthermore, by revisiting the config- 
urations discussed in the penultimate paragraph of the 
Introduction, we find that the kink energy £k becomes 
--i£i —eiiT- Thus, for a repulsive ear, |£i| will be larger 
than predicted by an analysis of, e.g., step-edge diffu- 
sivity that neglects ear- Lastly, the close-packed edge 
energy, i.e. the T = line tension /3(0) = — ^£1 —£2, be- 
comes — ^£1 —£2 —2€RT 



V. DISCUSSION AND CONCLUSIONS 

We now turn to experimental information about the 
interactions, followed by comments on the limited avail- 
able calculations of them, often recapitulating the dis- 
cussion in Ref. All the experiments are predicated on 
the belief that at 320K there is sufficient mobility to al- 
low equilibrium to be achieved. If the NNN interactions 
are to explain at least partially the high stiffness of ex- 
periment compared to Ising theory, the NNN interaction 
must be attractive and a substantial fraction of £1. Since 
compact islands do form on the Cu(OOl) surface, it is 
obvious that £1 is attractive. If £2 is also attractive, as 
required for reduction of the overestimate of ksT/ (3, then 
the low-temperature equilibrium shape has clipped cor- 
ners (octagonal-like, with sides of alternating lengths), as 
noted in Ref. |3 no evidence of such behavior has been 
seen. The lack of evidence of a decreasing stiffness near 
9 K, 45° suggests that £2/£i is at most 1/5. 

There is implicit experimental information for £2: from 
island shapes^'' and fluctuations^^ /3(0) = 220±llmeV. 
Since related measurements showed 3£i = — 128meV, we 
deduce £2 = —92 meV if eB.T is insignificant. These 
values imply that £2/£i is somewhat larger than 1/3, 
which seems unlikely in light of the unobserved predic- 
tions about the shape of islands in that case (cf. the end 
of Section III). 

To corroborate this picture, one should estimate the 
values of £1 and £2, as well as euT, from first-principles 
total-energy calculations. In contrast to Cu(lll))22i2^ 
however, no such information even for £1 has been pub- 
lished for Cu(OOl); there are, however, several semiempir- 
ical calculations which found £k ~ 0.14eV.^^ In such cal- 
culations based on the embedded atom method (EAM), 
which work best for late transition and noble fee met- 
als, the indirect ("through-substrate") interactions are 
expected to be strong only when the adatoms share com- 
mon substrate nearest neighbors; then the interaction 
should be repulsive and proportional to the number of 
shared substrate atoms. ■^^ (Longer range pair interac- 
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tions and multisite non-pairwise interactions are gener- 
ally very-to-negligibly small in such calculations; they 
probably underestimate the actual values of these inter- 
actions since there is no Fermi surface in this picture, and 
it is the Fermi wavevector that dominates long-range in- 
teractions.) If the NN and NNN interactions on Cu(OOl) 
were purely indirect, we would then predict €2 =561 > 0. 
However, whenever direct interactions (due to covalent 
effects between the nearby adatoms) are important, they 
overwhelm the indirect interaction. At NN separation, 
which is the bulk NN spacing, direct interactions must 
be significant, explaining why ei can be attractive. It is 
not obvious from such general arguments whether there 
are significant direct interactions between Cu adatoms 
at NNN separations. (For Pt atoms on Pt(lOO), the only 
homoepitaxial case in which ti was computed semiem- 
pirically, EAM calculations^^ gave e2/|ei| =0.2, less than 
half the ratio predicted by counting substrate neighbors, 
but with the predicted repulsive 62.) It is also not obvious 
a 'priori whether multi-atom interactions also contribute 
significantly. (For homoepitaxy, the only semiempirical 
result is that they are insignificant for Ag on Ag(OOl);^^ 
however, it is likely that semiempirical calculations will 
underestimate multi-atom interactions.) 

To address these questions, we are currently carrying 
out calculations^^ using the VASP packagei^Si Prelimi- 
nary results for Cu(OOl) suggest that £2 is indeed at- 
tractive, and that e2/ei is about 1/8; however, there are 
indications of a repulsive right-triangle trio interaction 
CRT with sizable magnitude (perhaps comparable to |e2|, 
consistent with a priori expectationsS^i^) , which would 
diminish rather than enhance the effect of 62- 

In summary, NNN interactions may well account for a 
significant fraction, perhaps even a majority, of the dis- 
crepancy between NN Ising model calculations and ex- 
perimental measurements of the orientation dependence 
of the reduced stiffness;^ the effect is even somewhat 
greater than estimated by the Twente groupSslSi. How- 
ever, inclusion of £2 is not the whole answer, nor, seem- 
ingly, is consideration of crt- One possible missing in- 
gredient is other multi-site interactions, most notably the 
linear trio eLT consisting of 3 colinear atoms (a pair of 
NN legs and an apex angle of 180°). In a model cal- 
culation their energy was comparable to eaxr^^ albeit 
with half as many occurrences per atom in the monolayer 
phase. The corrections due to e^T would be more com- 
plicated than simple shifts in the effective values of ei 
and 62- Since direct interactions are probably important, 
there is no way to escape doing a first-principles compu- 
tation; we continue to use the VASP package to extend 
our preliminary calculations.^* A more daunting (at least 
for lattice-gas afhcionadoes) possibility is that long-range 
intrastep elastic effects may be important. Ciobanu and 
Shenoy have made noteworthy progress in understand- 
ing how this interaction contributes to the orientation 
dependence of noble-metal stepsi^ 



APPENDIX: CALCULATIONAL DETAILS 

1. Partition Function 

To carry out the sum in Eq. we consider the 

Fourier transform of Z{Y): 



W(p) 



^exp^(f^Aj -if (A,)) 
{A} j=i 

00 

exp{ifiA-K{A)) 

,A=-oo 



, (A.l) 



where K {A) = {V + H \A\ + U \A - 1\ + D \A + 1|) is 
the energy in Eq. (^), associated with adjacent columns 
with height difference A. Carrying out the summation 
in Eq. (|A.ip gives 

= -^lnW{tfi) = V + U + D -\nB{tfi), (A.2) 



where 
B{in) = 1 



o2D 



^H+U+D+ifi _ I ' pH+U+D-ifj, _ I' 

Thus, the original partition function ZiY) is: 



(A.3) 



1 1"°° 

ZiY) = — df^e-^^'-Wifi) 



(A.4) 



1 

2^ 



d/i exp 



L —iiitanO — 



g(^A^) 



For i 3> 1, we can evaluate this inverse transform by 
steepest decent approximation. The saddle point occurs 
on the imaginary axis {fj, = —ip), at the value po given 
by the stationary-phase condition: 



g'(po) 
ksT 



= m — tan ( 



(A.5) 



Calculating the derivative from Eqs. ljA.2|) and (jA.Sp . we 
find 



m = B'{po)/B{po), 



(A.6) 



where prime stands for dp. The leading contribution to 
this integral (|A.4|I is just the integrand evaluated at this 
point: 



Z[Y) w exp 



L mpo + 



gjpo) 



(A.7) 
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2. Analysis of g"{p) and specialization to U — D 
From Eqs. HA.2p . we find 

g'ip) 



ksT 



'B\p)/B{p) 



(A., 



and 



^" = -B"{p)/B{p) + [B'{p)/B{p)f . (A.9) 



fcRT 



This can be simplified, by Eq. (|A.6|I . to 

9" iPo) 



= -mB"{po)/B'{po)+m^ 



(A.IO) 



the quantity needed for computing the stiffness as a func- 
tion of m. While straightforward, computing the deriva- 
tives with the general form for B (Eq. (|A.3() with p — ip) 
is quite tedious. A slight simplification emerges if we spe- 
cialize to the physically relevant case U = D. Then, with 
S = H + 2D,we have 



r,2D 



r,2D 



B{p) = 1 



= 1-e 



gS+p _ I e^^P — 1 
sinhS* 



■2D 



cosh S — cosh p 



1 - e^^ + 



2^^ sinh^ 
C{S,p) ' 



(A.ll) 



so that 



= e^^sinhS* 



sinhp 



(A.12) 



and 



= 6^^ sinh^ 



cosh p 2 sinh p 



_C2(^,p) C3(5,p) 

Inserting these expressions into Eq. (|A.6|I . we have 
sinh po sinh S 



(A.13) 



C(5, po) [sinh S - po) (1 - ' 
Similarly, with Eq. (|A.10p . we find 



9" jpo) 
ksT 



2 sinh Po 
C{S,po) 



coth Po 



(A.14) 



(A.15) 
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